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Abstract 

We study the Hamiltonian motion of an ensemble of unconfined clas- 
sical particles driven by an external field F through a translationally- 
invariant, thermal array of monochromatic Einstein oscillators. The sys- 
tem does not sustain a stationary state, because the oscillators cannot 
effectively absorb the energy of high speed particles. We nonetheless show 
that the system has at all positive temperatures a well-defined low-field 
mobility /i over macroscopic time scales of order exp(— c/_F). The mobility 
is independent of F at low fields, and related to the zero-field diffusion 
constant D through the Einstein relation. The system therefore exhibits 
normal transport even though the bath obviously has a discrete frequency 
spectrum (it is simply monochromatic) and is therefore highly non-Ohmic. 
Such features are usually associated with anomalous transport properties. 

Keywords: normal transport, inelastic Lorentz gas, diffusion, mobility 

1 Introduction 



The study of classical and quantum dynamical systems of a few degrees of 
freedom coupled to an environment containing many continues to attract con- 
siderable attention for various reasons. First of all, such studies shed light 
on fundamental issues of non-equilibrium statistical mechanics, such as return 
to equilibrium, the emergence of irreversibility and the existence of stationary 
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non-equilibrium states. More generally, they arise in the rigorous derivation 
of macroscopic transport laws (such as Ohm's law and Fourier's law) and of 
reduced equations of motion (e.g., Langevin equations and master equations) 
from a microscopic, Hamiltonian description of the dynamics. The explicit com- 
putation of transport coefficients such as the diffusion constant or the mobility 
from first principles is another motivation of such studies, particularly in solid 
state physics. It is this problem that interests us here. 

The case in which the environmental variables to which the particle couples 
are vibrational degrees of freedom is of obvious interest to many fields of physics 
and is the most amenable to both analytical and numerical studies. The model 
Hamiltonian we will study in this paper was introduced in [24j and can be 
written 

Here q and p are the particle position and momentum and {qrmPm) are the 
position and momentum of an oscillator of frequency w located at 
with which the particle interacts provided \q — ma\ < a; indeed, p vanishes 
when \q — mal > a and equals unity otherwise (we take a ~ 2a = L > 0) . Here 
> is a driving field and we are interested in computing the low field mobility 
of the particle as a function of the temperature of the oscillator bath. 

Our motivation for studying this particular model is twofold. On the one 
hand, the model obviously resembles those describing electron-phonon interac- 
tions in solids [171 El] , an analogy that is further worked out in [11] . Alterna- 
tively it can be viewed as describing the motion of a particle in an inelastic one- 
dimensional Lorentz gas. Recall that in the usual Lorentz gas, a popular model 
for the study of transport properties, a particle moves on a two-dimensional 
plane on which circular hard scatterers are randomly or periodically placed. In 
order to describe the interaction with a thermal environment, the scatterers can 
be allowed to rotate, as in [T3ll2Q], yielding a model similar to ours but with the 
environment described through rotational degrees of freedom rather than vibra- 
tional ones. Alternatively, the scatterers can be kept fixed, as in [5117], while the 
environmental interaction is phenomenologically described through the use of a 
Gaussian thermostat; in the latter situation only the particle degrees of freedom 
remain dynamical and the system is no longer Hamiltonian. Our model can be 
viewed as a one-dimensional Hamiltonian version of such models. Unlike previ- 
ously studied two-dimensional transport models, however, the present model is 
not homogeneous in energy; as a result it displays a wide range of microscopic 
behavior that varies with the temperature of the system. 

The main conclusion of the paper regarding the present model is that, pro- 
vided the driving field F is low enough, the particle distribution quickly acquires 
a drift velocity vp proportional to the driving field F which it then maintains 
over times that are exponentially long in F^^. This current is metastable, how- 
ever, due to a runaway phenomenon that arises because oscillator lattices, in 
general, do not efficiently slow down very fast particles (see Section [5]). Thus, 
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this system does not approach a stationary non-eqiiihbrium state. We nev- 
ertheless show that the associated mobility is well-defined and related by the 
fluctuation-dissipation theorem to the diffusion constant D and the inverse tem- 
perature [3 = 1/kT through Einstein's relation /i = (3D. Thus, in this simple 
fully Hamiltonian model, the particle exhibits normal transport behavior obey- 
ing Ohm's law at low fields. 

These conclusions will be surprising to some, since previous studies of sim- 
ilar models seem to indicate rather strongly that coupling to a monochromatic 
bath would be expected to lead to anomalous particle transport. To support 
this claim we conclude this introduction with a short overview of the relevant 
literature. There is a considerable body of work on a variety of models in which 
a particle of momentum p and position q couples linearly to the displacements 
or momenta of a heat bath of oscillators. A very general class of such models 
was proposed in ^6]. To describe the results in the literature relevant to our 
work here it will be sufficient to consider a specific subclass of models described 
by the following Hamiltonian [5]: 

+ [ [ p(x - q)c{uj(£,))q(x,0<yi{x)<J2{£.)dxd^. (2) 
J D J B 

Here X — (?(a;, i^),p(a;, ^))2,g£)^K,4e-BcR'* where q{x,£,), p{x,^) are the displace- 
ment and momentum of an oscillator of frequency uj{^), which should be thought 
of as finding itself at the point x € D. The parameter ^ is to be thought of 
solely as an index for the oscillators. As a result, the collection of oscillators 
(g(a;,^),p(x,^))^gBcR''j ^ being fixed, represent a bath of vibrational degrees of 
freedom of the medium at the point x € D. Furthermore cri(a;)dx and (T2(^)d^ 
are measures (on D and B respectively) and p{x), c(w(^)) positive functions de- 
termining the coupling of the oscillators to the particle. Note that, since neither 
(72 nor c depend on x, the oscillator baths have the same characteristics at all 
points X ^ D. We will choose for p a compactly supported positive function 
centered at the origin. The Hamiltonian ([1]) is of the above type, with 

^^S{x - 'ma), (72(0 = '^l'? - w), 

m 

and uj{^) = uj. 

Previously studied models with Hamiltonians as in ^ are of two types. 
Models with a continuous distribution of oscillators uniformly throughout all of 
space are studied in [H 12] ; in that case ai (x) = 1 and D = R. Alternatively, in 
[HIH], the oscillators are located at a countable set of points D = {xm\m G Z} 
chosen randomly and homogeneously throughout space so that 

(Jlix) ^ ^S{X- Xm)- 
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To discuss the results obtained in those previous studies we introduce the 
couphng weighted spectral density of the system, defined as 

7rc(w)2 

where the spectral density is defined via 

JlAj{^)<UJ Jo 

In [U [21 m [9] , the distribution of frequencies is supposed to be continuous and 
~ w** as ^ 0, for some s > 0. It is then argued in [U HI |9| that 
normal transport properties are obtained in such models only if s = 1. In 
particular, these authors argue that when V{q) = 0, particle motion is diffusive 
provided s = 1, whereas anomalous diffusion occurs otherwise. Similar results 
are obtained rigorously in p] where it is proven that in the fully translationally 
invariant models considered (namely cri(x) = 1) a particle moving at constant 
speed V through an oscillator bath at zero temperature experiences a friction 
force proportional to at low v: this again indicates that a finite mobility 
will be obtained only if s — 1, a fact rigorously proven in pi for the model 
considered, at zero temperature. 

Similar conclusions have been obtained in the much more numerous studies 
([21 HI 13 [El El [H [2l] and references therein) dealing with a simpler but 
related class of models in which the potential V (q) is chosen to be confining. 
Note that this prevents the particle from making arbitrarily large excursions 
from equilibrium, justifying a "dipole approximation" in the Hamiltonian ^ 
in which the coupling is linear in q as well. After some rewriting, the resulting 
bilinearly coupled Hamiltonian takes the form 

HcL{q,P,X) = P^ + Viq) + ^l^{p{0'+u;iO'qma2m^ 

+q [ c(t^(C))g(C)^2(0de 
Jb 

and is then much simpler to analyze. In classical systems, for example, it 
straightforwardly leads to a generalized Langevin equation for the particle: 

q{t)+ I l{t-s)q{s)ds = -^V,s{q{t)) + f{t), 

with 

VMl) = Viq) - \ 

and where 

Jb ^ (i) Jq ^ 



Jb ^ {Q 
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is the retarded friction kernel and 



fit) = - / iu^iO) 



qoiO + ^qo) ] cos(a;(C)t) + MO sin(a;(C)i) 



(3) 

the random force exerted by the medium on the particle, which depends on 
the initial conditions for the particle (9(0) = q^) and for the bath variables 
(9o(C)iPo(?))- Return to equihbrium for such systems was shown for general 
confining potentials in [18] provided s < 1; the particular case s = 1, with a 
harmonic potential V was treated in O [15] . Note furthermore that the same 
model, with V{q) = can also be understood as a model for a non-confined 
particle: in that case the oscillator degrees of freedom can be thought of as 
internal degrees of freedom of a composite particle, rather than as environmen- 
tal degrees of freedom. For this model, it was shown in [T21 UHl US] that the 
particle diffuses only in the Ohmic case s = 1 and is superdiffusive (respectively 
subdiffusive) when s > 1 (respectively s < 1). 

In short, in all existing literature in which theoretical results are available 
on the behavior of a particle coupled (linearly or not) to vibrational degrees 
of freedom, normal transport properties for the particle occur only if the os- 
cillator bath has a continuous oscillator spectrum with a sufficient density of 
low frequency modes (s = 1). Our study here shows that while an Ohmic local 
bath is certainly a sufhcient condition for normal transport, it is not a necessary 
one. On the contrary, we show that even with a monochromatic oscillator bath 
having no low frequency modes, a finite diffusion constant and well behaved low 
field mobility can be obtained. 

The paper is organized as follows. In Section [2] we describe our numerical 
results on the particle mobility and its temperature dependence. Section [3] 
provides an explanation of the high temperature part of our data in terms of 
a random walk model motivated by the underlying physics of the Hamiltonian 
dynamics. In Section [5] we apply linear response theory to prove an Einstein 
relation for the model, valid at all temperatures, which we apply to explain 
the low temperature behaviour of the model. In Section [5] we show that no 
stationary state is established in the system due to a runaway phenomenon 
inherent in all models of this type. 



2 Mobility: numerical results 

A simple inspection of the equations of motion shows that, under the dynamics 
generated by the Hamiltonian ([1]) , the particle undergoes at all times a constant 
acceleration F. Its motion is decoupled from the motion of the oscillators, except 
at those instants of time when it reaches the edge of one of the interaction 
regions, in other words when \q{t) — ma\ = a for some to e Z. At those times, it 
encounters a potential energy step A of magnitude |a(7„i| and it therefore either 
turns back or proceeds with a change in momentum that is simply determined by 
energy conservation. If the kinetic energy of the particle at the moment it arrives 
at the edge of the interaction region is greater than the potential discontinuity. 
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Figure 1: Mean square displacement of the undriven particle {F = 0) as a 
function of ujt for several sets of parameters and temperatures. In both panel (a) 
(where ci = 0.015, C3 — 0.5) and panel (c) (where ci — 0.5, C3 = 0.5) triangles 
indicate C2 — 0.5 and circles C2 — 5. In both panel (b) (where ci = 0.002, 
C3 = 2,) and panel (d) (where ci = 0.07, C3 = 2,), diamonds indicate C2 = 0.5 
and squares C2 — 5. Each curve is an average of 10'^ — 10^ trajectories. Straight 
lines are linear fits and ci, C2, C3 are defined in (H])-®. 

it passes into the next region with a new kinetic energy reduced or increased by 
the potential step that it traverses in doing so. Otherwise it undergoes an elastic 
reflection. The oscillator at Xm = ma, on the other hand, oscillates about the 
displaced equilibrium = — a/w^ at those times t for which \q{t) — ma| < a, 
and about the unperturbed equilibrium = otherwise. 

These features allow for an event-driven numerical implementation of the 
dynamics, without the need for numerical integration of a non-linear differen- 
tial equation. As a result, we have been able to run simulations over very long 
times (of the order of 10^ — 10^ oscillator periods), computing average displace- 
ments or mean square displacements for the particle, with an ensemble of initial 
conditions for the particle-oscillator system drawn from the (field-free) Boltz- 
mann distribution. To obtain satisfactory statistics, we ran between 10"^ and 
10^ trajectories for each set of parameters considered. 

In this paper we are concerned both with the high and the low temperature 
regime. To understand what this means, note that there are two mechanical 
energy scales associated with this model. The first one is the binding energy 
—Eb = — a^/2a;^, which is the classical ground state energy of the system, 
obtained when the particle is at rest in one interaction region (for example the 
one associated with the oscillator at x = 0), with the corresponding oscillator 
in its displaced equilibrium position and all other oscillators at rest (q™ = = 
Pm,m 0). The other one is the kinetic energy £'0 — a'^LiP of a particle that 
crosses the interaction region of length 2a in one oscillator period. Along with 
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Figure 2: Mean particle displacement as a function of ojt for several sets of 
parameters, temperatures, and forces {Fq = auP'). In both panel (a) (where 
ci = 0.015, C3 = 0.5) and panel (c) (where ci = 0.5, C3 = 0.5) triangles indicate 
C2 = 0.5 and circles C2 = 5. In both panel (b) (where ci = 0.002, C3 = 2,) and 
panel (d) (where ci = 0.07, C3 = 2,), diamonds indicate C2 = 0.5 and squares 
C2 — 5. Each curve is an average of 10'^ — 10^ trajectories. Straight lines are 
linear fits. 



the thermal energy kT — 1//3, this allows the introduction of two energy related 
dimensionless parameters, which we denote by 

ci EE (3Eb and C2 = Eb/Eq. (4) 

The model also has two length scales L ^ a — 2a and cr, leading to a third 
dimensionless parameter 

C3 = 2a /L. (5) 

By high temperatures we will mean the regime ci << 1, (3Eo = C1/C2 << 1, for 
which the average thermal energy of the particle is (much) higher than both the 
binding energy Eb and Eq. In the low temperature regime, these inequalities are 
reversed. As we will show below, the temperature dependence of the mobility 
is in this model qualitatively different in these two different regimes, due to the 
qualitatively different nature of the microscopic dynamics. 

In Fig. [T] plots are seen of the numerically computed mean square displace- 
ments of the particle for various parameter values, in absence of a driving field. 
As already observed in [21], the motion is diffusive. 

In Fig. [5] we present plots of numerically calculated values of the averaged 
particle displacement {q (t)) in units of a, in the presence of a driving field for 
several systems and temperatures. In each case, the average is computed over 
initial conditions of the particle-chain system which are drawn from its joint 
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Figure 3: Drift velocity vp as a fiinction of F. Straight lines are linear fits. 
For all, C2 — h, vq = ato and Fq = auP' . In the top (bottom) panel, triangles 
indicate ci = 0.5 (0.015), C3 = 0.5; squares indicate ci = 0.9 (0.03), C3 = 1; 
circles indicate ci = 0.7 (0.02), C3 — 2. 



thermal equilibrium distribution in absence of the driving field. One observes 
that {q [t)) is linear in t over macroscopically long times and distances. Over 
the time scales explored, each system investigated thus attains a well-defined 
drift speed vp, that is easily determined through linear regression. 

As seen in Fig. |31 for fixed T and system parameters the resulting drift 
speeds are, for sufhciently small fields, linear in F. This allows us to extract 
through linear regression, a well-defined low field mobility which we have plotted 
in Fig. [4] as a function of (SEb ■ 

At low temperatures, the numerics is complicated by the following phe- 
nomenon. First, since we compute displacements and since the system is peri- 
odic, we can start all particles in the cell &i q — Q. Then 
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Figure 4: Numerically computed mobilities n/n^ (see ([9]) for /i^) as a function 
of (3Eb over a wide range of temperatures and for nine different sets of model 
parameters, distinguished with different symbols. For each parameter set, the 
mobility is shown for six different temperatures. The straight line corresponds 
to ri {PEbV^''^ with 7? = 0.32. 



are the probabilities that the particle finds itself respectively in the interaction 
or in the non-interaction region of the cell [2l]. Note that, obviously, the latter 
is exponentially less likely as the temperature is lowered. In addition, only 
an exponentially small fraction of the particles in the interaction region will 
ever leave the first cell under the dynamics, because most particles will never 
manage to overcome the typical energy barriers they meet at the edges. This 
leads to an exponential suppression of the mobility at low temperatures that 
can be observed both on Fig. [4] and Fig. [5] and that will be further explained 
in Section 21 These phenomena also affect the numerics in the following way: 
trapped particles will of course not contribute to the mobility and running the 
corresponding trajectories would be a waste of computer time. It would also 
keep one from getting good enough statistics, especially at low temperatures, 
where the drift speed is low. Fortunately, using the analysis of lOj, it is possible 
to predict a priori from the initial conditions of the particle and the oscillator if 
the particle will ever get out, and implementing this in the program allows one 
to avoid running the corresponding trajectories uselessly. It follows from this 
analysis that particles can remain trapped only if they start in an interaction 
region. Suppose therefore {q,p) are the initial data for the particle, with < 
g < 2f7 so that the particle is indeed in the interaction region, and let (qojPo) 
be the initial conditions for the oscillator in the cell at the origin. The total 
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energy of this particle-oscillator system can then be written 
E = e(g) + £080(90,^0) - Eb - 2Fa, 

where 

eosc('7o,Po) = ^(Po + (90 + of)), £{q) = - F{q - 2cr). 

Analyzing the particle-oscillator dynamics along the lines of [10] the necessary 
and sufhcient conditions guaranteeing that the particle will leave the cell at the 
origin at some point in the future are 

2Eb ~ e{q) > eosc(go,Po) > j^i^il) ~ '^Esf . 

We now turn in the next sections to the theoretical explanation of the ob- 
served numerical data described here and in particular of the temperature de- 
pendence of the computed mobilities observed in Fig. 2] and Fig. [51 

3 High temperature mobility: a biased random 
walk model 

In order to explain theoretically the high temperature part of the numerical 
results described in the previous section, we will further develop the simplified 
random walk model of the fully Hamiltonian dynamics at high temperature 
that was developed in [24] for the case F = 0. This will provide us with a 
clear and simple physical picture of the motion of the particle in this regime. 
First note that, at high T, the typical kinetic energy of a particle (which is of 
order kT) is much higher than the typical barriers it encounters (which are of 
order y/EBkT). As a result, particles drawn from a thermal distribution tend 
to enter most interaction regions they encounter, rather than reflecting. Now 
any particle with a high velocity v will cross the interaction region in a time 
a /v which is short compared to the oscillator period and as a result it will tend 
to lose an energy AE = —AEbEo/v'^ to the oscillator, as shown in [24 . At 
high temperatures, one has PEq << 1 and a particle with approximate thermal 
speed will satisfy the above condition since vth — l3~^^^ so that a/vth « l/uj 
is equivalent to Eo(3 « 1. 

On the other hand, the same particle gains or loses an energy Fa from the 
field while traversing the cell, depending on whether it is moving to the right 
or to the left. Consequently, if the field is too large. Fa > AEbEq/v^^, particles 
with thermal speed moving to the right will be accelerated indefinitely, so that 
it is not possible for a drift speed to be set up in the system. At small fields, on 
the other hand 

Fa« 4Eb Eo/vl^ = A/SEb Eo (6) 
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such particles will lose a net energy. The resulting average velocity change upon 
traversal of one cell is given by 



Av = ApE- 



.1 



AEbEq 



1 



-Fa. 



(7) 



So, when the temperature is high and the field small, the energy loss experienced 
by most particles in the thermal distribution is small, and they traverse many 
cells before slowing down to a point where their kinetic energy is smaller than or 
comparable to the typical barrier heights. At that point, they can be expected 
to receive a randomizing "kick" from one of the oscillators. We will assume 
that this randomizing "kick" reinstates a thermal distribution for the particles. 
One easily finds the distance £f{v) over which a particle of high initial velocity 
V travels before slowing down, as well as the characteristic time Tp{v) of the 
slowing down process, by integrating ([7]): 

7l3 1 

r-r dv, 



and 



tf{v) 



— av — 
dv 



4:EbEo 1 - 



1 



4Eb Eq 



V\V\ 



V AEbEq 1 



Fa . 

iEgEo 



dv. 



This yields: 



and 



Tp{v) = t(v) 



£f{v) =e{v) 



3 Fav^ 



1 



5 iEgEo \v\ 



Fav^ 



t{v) 



12EbEq 



e{v) = -VT (v) . 



GEbEo \v\ 

Note that, as expected, for v > 0, one has £f{v) > £f{—v)- particles moving to 
the right are pushed on by the field and therefore travel further before slowing 
down. 

The result of this analysis is that at high temperatures one can view the 
particle motion as a biased random walk with steps £f{v) and waiting times 
tf(u), with V drawn after each randomizing kick from the Boltzmann distri- 
bution (/?/27r)^/^ exp(— /?'i;^/2). This picture therefore leads one to predict an 
average drift velocity 



Vf 



{If 



Fa 



{tf) SEbEo 



and thus a mobility 




-3/2 _ 



-3/2 



(8) 



(9) 



The data presented in Fig. |4] show that this simple picture of the dynamics 
yields the correct temperature behaviour of the mobility at high temperatures 
as well as a reasonable estimate of the prefactor. 
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Figure 5: A log- linear plot of fi/^J-l (see (fT2|) for the definition of as a 
function of PEb for the same model parameters as in Fig. [4l emphasizing the 
low temperature data. The curve corresponds to 0.14(2:^/^ + l.5x~^'^^^)exp{—x) 
with X = (3Eb- This clearly brings out both the exponential dependence at low 
temperatures and the power law dependence at high temperatures. 



4 Low temperature mobility: linear response 

It was shown in [24j that the physics underlying the low-temperature dynamics 
of the model at zero field {F = 0) is of a very different nature than the one 
at high temperatures. It can again be described in terms of a random walk 
but this time with nearest-neighbour hopping of the particle. The analysis of 
[24j yields a very good prediction for the diffusion constant and its temperature 
dependence at low T: 



D ^ D 



exp {-/3Eb) 



th" 



where 



*^ 2a \ 27r 



f3EB > 1 



(10) 



(11) 



In particular, it is seen from this result that the diffusion constant is thermally 
activated due to the self-trapping of the particle at low temperatures briefly 
explained in Section [5] (see ^24j for details) . This low temperature analysis is 
however more involved than the high-temperature one and will not be further 
adapted to the non-zero field situation we consider here. Instead we will present 
in this section an alternative derivation of the mobility of the system and of its 
temperature dependence via Kubo's linear response theory. We will show this 
leads to a finite time variant of the Einstein relation valid at all temperatures 
and for small non-zero fields which, combined with (|10p . yields for the low 
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temperature behaviour of /i the result 



^Je^ exp {-PEb) , Ml = ^ V^E^ ^^^^ 

This is in very good agreement with the low temperature numerical data as can 
be seen in Fig. [51 

Let us consider again the Hamiltonian in JT]) and write it as 

HF{q.p-,X)^Ho{q,p-X)-Fq, 

where X — {qm.,Pm)m£Z- We will write $f for the phase space flow generated 
by this Hamiltonian and use the shorthand {q,p; X) = {q{t; F) , p{t; F)) . In 
other words, q{t, F) and p{t; F) are the position and the momentum at time t 
of a particle that at time had position q and momentum p and that evolved 
under the dynamics generated by Hp in a medium with initial condition X. 
We will be interested in the thermally averaged mean velocity of the particle at 
time t, 

{iqit;F)-q)) 
VF[t) = , 

which is the thermal average of its mean displacement, divided by the time t. 
Also, we define the mobility by 

Note that this is, a priori, a time and field dependent quantity and that vpit) 
and fJ^pit) are precisely the quantities that we compute numerically. We are 
interested in their long time and small field behaviour. Here and in what follows, 
for any phase space function f{q,p; X), (/) designates the thermal average of / 
for the Hamiltonian Hq. In addition 

A/t = (/o$f >-(/>. 

A formal application of first order perturbation theory to our system then shows 
readily that for any phase space function /, 

Aft f ds F({g,iJo}/o + Ot{F^). 

Jo 

Consequently, if / = g, one finds, after a change of variables in the integral 

VF{t)^PF- [ ds {pq{s;Q)) + Ot{F^). 
t Jo 

Writing furthermore 

(j(s;0) -q= ds'p(s';0) 
Jo 
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Figure 6: Plot of i.i/{(3D) against (3Eb for the nine systems used in the previous 
figures. Six temperatures are represented for each system. 



this can be rewritten, after changing the order of the time integrations: 







fds' 




'o 





On the other hand, let us introduce the mean square displacement of a freely 
evolved particle in thermal equilibrium 



m := 



{{q{t;0)-qy 
2t 



This can be rewritten in the usual way as 

D{t)^ j\l~'-){pp{s)) ds. 

Hence, we immediately obtain the following finite time version of the Einstein 
relation: 

Mi)^PD{t) + Ot{F), (13) 

which holds for all times and fields. 

For fixed time the error term is order F , but in view of the metastability 
of the current alluded to in the introduction and further explored in the next 
section, one cannot expect this term to remain small uniformly in t. In other 
words, the limit as t goes to infinity cannot be taken in (fT^ . In fact, the system 
does not sustain a stationary state. On the other hand, the results of [H] (see 
also Fig. [T|) strongly indicate that 



D := lim Dit) 



(14) 
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Figure 7: Log-log plot of {q{t,F))/F against t for several values of F/Fq, as 
indicated and for the parameter values i3Eb = 0.5, Eb / Eg = 5,2a /L = 0.5. 
The parts of the graphs parallel to the dashed line correspond to behaviour 
linear in time (slope equal to 1). The common j/-intercept with the vertical axis 
reflects their common value for the mobility. But for all forces, at large enough 
times, the displacements increase faster than linearly with time. 



does exist, for all temperatures. In fact, it is clear from the numerics that the 
diffusive regime is reached on a short time scale. At such times, it suffices to 
take F small enough to make the error term in (|13p small, leading finally to the 
following finite time version of the Einstein relation, 

^ipit) ^ f3D{t) ^ pD, 

linking the numerically computed diffusion constants and mobilities. Note that, 
since the right hand side depends neither on t, nor on F, this means that, for 
such times and fields, the same must be true for urit)- We tested this relation 
numerically for nine systems, each at six different temperatures, as shown in 
Fig. El 



5 Runaway 

As we have already suggested, the system studied here cannot sustain a constant 
drift speed vp indefinitely. This can be observed by considering quite strong 
fields, as we have done in Fig. [71 in which one clearly sees the drift speed 
increasing rather suddenly at some critical time tc{F) that depends on F. We 
will analyse this phenomenon below in order to obtain an estimate for tc{F). 
It is clear from dZ]) that, on average, particles with velocity 

V > y^4:EBEo/Fa 
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gain more field energy than they dissipate, so that their mean velocity increases 
eventually as v (t) ?' + Ft. In addition, one may expect that even somewhat 
slower particles ma.\. us a result of thermal fluctuations, see their speed in- 
crease beyond ^/4:EBEo/Fa and then accelerate indefinitely. This suggests the 
definition 

Vc{F) = ^^AEeEo/Fa (15) 

of a critical speed Vc{F) where 7 is a numerical constant of order 1. Particles 

with speeds larger than Vc{F) will accelerate indefinitely, on average. Note that, 
obviously, this critical speed increases when the applied force decreases, and 
is otherwise independent of the temperature. Now, in a thermal distribution 
there are always particles with speeds larger than Vc{F) and their increasing 
velocities will eventually dominate the average velocity, which will no longer 
remain constant. Provided the field F is low enough, however, the thermal speed 
Vth = (3~^^'^ is much lower than the critical speed Vc{F). Particles destined to 
undergo runaway are then rare, since they find themselves in the tail end of the 
distribution. Indeed, 

Vib < v,,{F) 4^ Fa < h^EeEa (16) 

and hence the fraction of particles in the thermal distribution that have a ve- 
locity V > Vc{F), which is given by 

PF = (^) ' r dv ~ {8nP^^EBEo/Fa)~^^\xp {-2I3^''EbEq/ Fa) 

is exponentially small in Consequently, for low enough fields most particles 

have speeds less than Vc{F) and to them the arguments of the previous sections 

apply, on average: they reach a limiting drift velocity vp- 

The average velocity of the ensemble will, therefore, take the form 

{v (t)) ~ vf{1- pf) + \U- {P + Ft) e-^f'/2 dp 
~ vf{1- Pf) + — i=e-5^'=(^)''^ + FtpF 



g 2 2EbE0 

Vf + [Fa) ' 



Consequently, for times t long enough for the slow particles to reach their 
final average velocity, but short enough that the second term in this last relation 
is negligible, the particles will exhibit a finite and constant current, and thus 
an Ohmic response. Thus, defining a critical time tc{F) when the fast particles 
start to dominate the mean particle current [i.e., when pFFtc{F) t;ir], we find 
that 

to{F) = VFa^ ^3^3 e (17) 
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Figure 8: Log- log plot of the critical breakdown times as a function of F for 
the system with parameters f3EB = 0.5, Eb/Eq = 5, 2a/ L = 0.5. The circles 
on this graph indicate numerically established values of (F) . The full curve is 
tc{F) where 7 — 0.6 and as in the previous figures Fq = a^w. 



which is exponentially long in F^ . 

We note that such an effect will be unobservable if the length of the sample 
supporting the current is much shorter than the exponentially large distance 
Lp = vptciF) the particles would drift during this time in an infinite sam- 
ple. We have tested this runaway phenomenon numerically for the system with 
PEb = 0.5, Eb/Eq — 5, 2a/L = 0.5. In Fig. [7] one observes for this system 
the numerically computed value of {q{t))/F for five values of the force F/Fq 
between 1 x 10"'^ and 2.5 x 10"'^. One observes that the displacement grows 
faster than linearly beyond a breakdown time, indicated with vertical arrows on 
the figure. This breakdown time is of the order lO'^w"^ — lO^w^^ for the three 
largest forces used and is seen to increase with decreasing force. The numer- 
ically determined values of this breakdown time are plotted against the force 
F/Fq, in Fig. [5] One notes that the breakdown time increases very fast with 
decreasing fields, in a manner qualitatively compatible with the exponential law 
derived above. Indeed, the full line in Fig. [8] corresponds to tc{F) in p7|) with 
7 = 0.6. This is compatible also with the following observation: for much lower 
forces, of the order of lO^^^i^o - 10~^Fo, the mean displacement of the particle 
increases linearly in time for the same system as seen in the circular data points 
of panel (c) in Fig. [2] for times up to i = lO^w"^. This strongly indicates that 
for those very small forces, the breakdown time will be considerably larger than 
lO^w-i. 

The previous considerations make clear what the difficulties and practical 
limits are of the numerical computations on this model, specifically at high 
temperatures. Indeed, because of the existence of the breakdown time, we have 
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to use very small external forces to study the mobility (see (HH])). But then the 
drift will be small and so the motion is diffusion-dominated for quite long times. 
Consequently, to observe the drift, one has to run many trajectories for long 
times in order to obtain sufficiently good statistics. 

The analysis of the runaway phenomenon allows for a closer study of the 
long time and low field limit of (jl3p . First, since this is a first order formula in 
F, at fixed t, and in view of (I14p . it is clear that 

u := lim lim upit) — 0D 

exists, with the two limits taken in that order. Due to the runaway phenomenon 
exchanging the limits is not possible: one cannot expect the limit as t — > +00 
of the left-hand side of to exist, for any fixed value of F, however small. In 
other words, this model does not sustain a stationary state. Nevertheless, our 
analysis also shows that the following precise mathematical statement can be 
expected to be true for this model: 

^^^t)^^^f3D + 0{F)+ei^) 
where the error term in F is uniformly small for all times in the range 

fi «t « tc{F). 
In other words, the following limit should exist: 

n = lim Upit) ~ 0D. 

M<<t<<tc(-F) 

6 Conclusion 

In conclusion, we have shown in this paper that an array of monochromatic 
oscillators at a positive temperature can serve as an effective heat bath for a 
particle driven through the lattice under the influence of an external driving 
field F, inducing for the particle normal transport properties over very long 
times and distances and this in spite of the fact that a truly stationary state 
cannot be sustained by the system. 

To end, we point out that the situation here is therefore different from the 
one proven to hold for the Lorentz gas with a Gaussian thermostat jSJ [Zj ■ In 
that case a stationary state is shown to exist and the limits in (jl3p can be 
exchanged. But there, friction is provided by a Gaussian thermostat and the 
model is therefore not fully Hamiltonian, as pointed out before. On the other 
hand, in the Lorentz gas with rotating scatterers [T31 [50] , the Einstein relation 
is checked numerically to hold, but for rather short chains and times. Finally, in 
[11] we have argued that replacing the kinetic energy of the particle by a 
cosine band Hamiltonian 2V^(cosp— 1) will suppress the runaway phenomenon, 
so that in that case as well a stationary state will exist. 
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